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' _^ ' Abstract 

String theory/M-theory generally predicts that axionic fields with a broad 
^ ^ mass spectrum extending below 10~^''eV are produced after compactification 

U to four dimensions. These axions/fields provoke a rich variety of cosmophysical 

phenomena on different scales depending on their masses and provide us new win- 
dows to probe the ultimate theory. In this article, after overviewing this axiverse 
idea[T], I take up the black hole instability as the most fascinating one among 
^ such axionic phenomena and explain its physical mechanism and astrophysical 

predictions. 

^ 1 Introduction 

00 

At present, supcrstring theory/M-theory appears to be most close to the ultimate 
T-H theory of Nature. Therefore, it is crucially important to find a clue indicating that 

these higher-dimensional theories are really behind our universe. 
^ Because the ultimate theory is a kind of UV completion of our low energy effective 

t ""j theory, its characteristic new features such as the existence of extra-dimensions in 

rS general show up in high energy phenomena. Thus, one natural approach to probe the 

ultimate theory is to study phenomena at the high energy end experimentally. Collider 
experiments like LHC and ILC are examples, but it is very difficult to raise the maximal 
energy to string scales in the near future. From this respect, it is more promising to 
probe inflation through CMB and gravitational waves, and look for cosmological relics 
from the early universe such as dark matter, gravitational wave background and cosmic 
strings, directly or indirectly. 

We can also probe the ultimate theory through low energy experiments. For ex- 
ample, if the size and structure of the extra dimensions is not completely stabilized, 
the values of fundamental constants may vary in time or spatially on cosmological 
scalesOinHUIS]. Further, there may exist light moduli fields. If they contains a scalar 



^This article is basen on the lecture given at the 2011 Shanghai Asia-Pacific School and Workshop 
on Gravitation, and will be published in its proceedings as a volume of the International Journal of 
Modern Physics Conference series. 
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field, it mediates a new force with a range corresponding to the Compton wavelength 
1//Lt. In fact, lots of experiments have been done to look for such a new force in submm 
ranges, although only upper bound have been obtained so far. [6] 

In contrast, if these light moduli fields are pseudo-scalar, i.e., axions, it is difficult 
to detect them by new force search experiments, because the force mediated by such a 
field is proportional to the spin or velocity of matter source and decreases faster than 
with the distance r from the source even when r is shorter than 1/Ai.[7] Hence, 
the mass of the axion can be very small. In particular, the Compton wavelength 
can become of the order of astrophysical objects or cosmological scales. In such a case, 
axion fields may provoke cosmophysical phenomena, as systematically discussed by 
Arvanitaki et al in the axiverse paper [[T]]. In this article, we take up the superradiant 
instability of black holes and astrophysical phenomena provoked by it as the most 
fascinating one among various new phenomena provoked by superlight axions. 

This article is organized as follows. First, in the next section, we briefly overview 
the axiverse idea and the cosmophysics based on it. Then, we focus on the black 
hole problem. After reviewing the basics on black holes, we discuss the superradi- 
ance by a rotating black hole, the superradiance instability of an axionic field and its 
astrophysical implications in order. 

2 String Axiverse 

As is well-known, axion was first introduced into physics as a pseudo-Goldstone boson 
for the Peccei-Quinn symmetry to resolve the strong CP problem[S]. This QCD axion 
was originally assumed to have a mass of MeV order and interact rather strongly with 
particles in the standard model. However, it was soon recognized that the existence of 
such a particle contradicts experiments unless its coupling to quark unless unless its 
coupling to the SM sector is extremely small, i.e., it is practically " invisible" . [3] 
The basic features of this invisible axion are summarized as follows. [TU] 

1. a neutral P- and CP-odd scalar coupled very weak to matter : 

gaqa{q-f5q) : gaq « niq/ fa, fa ^ lO^GeV, (la) 
QajaFAF: ga^-l/fa (lb) 

2. SmaU mass by the QCD instanton effect: - 10-3eV(10i°GeV//a) 

3. Dark matter candidate: ^a = 0.01(/a/10i°GeV)i-i^5. 

Here, note that the interactions ([T]) have shift symmetry, i.e. they are invariant under 
the transformation a — a -f const mod. field equations. This feature has a crucial 
importance in protecting the axion to get a large mass by various quantum corrections. 
On the basis of this, the concept of the axion is now generalized to include all pseudo- 
scalar particles/fields that have P and CP violating interactions with shift symmetry 
at the tree level. 
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2.1 String axions 
2.1.1 Origin 

The axiverse idea is based on the fact that axions defined above are expected to be 
produced abundantly by realistic compactifications of string theory /M-theorv|l 11 [T]. 
They come form fields contained in any string theory as essential ingredients|12j. 

For example, the bosonic sector of the heterotic SST contains, in addition to the 
spacetime metric g, a dilaton (p and gauge fields A, a 2-form potential B whose action 
in the string frame reads 

2kIoSb^ ~-e-''^*HAH; H ^ dB ^ ^ {u^s - ^hs) , (2) 

where ol is the inverse of the string tension, and w^g and w^g are the Chern-Simons 
connections for the gauge field A and the gravitational connection w, respectively: 

dwgs = Tr (FAF), duj^s = tr A ^) . (3) 

This action is invariant under the gauge-transformation 

M = dA, Sio^^dQ, SB ^da + {a'/A){TT{XdA)+tT{Qdio)} (4) 

with an arbitrary 1-form a. 

When the theory is compactified on a Calabi-Yau 3-fold F to a four-dimensional 
spacetime X, the form field B produces two types of axionic fields in X. To see this, let 
us consider a simple product-type compactification, ds^{Mio) = ds^{X4) + ds'^{Yia). 
Let 77* (i = I,-- - ,b2{Y)) be a basis of harmonic 2- forms on Y dual to a basis of 
H2{Y,Z). Then, B can be expanded as 

b2{Y) 

B = llY,adxX + P{x), (5) 

i=l 

where ts — Itt^/o' and /3(x) is a 2-form on X^. By inserting this to the original action 
S'b, we obtain 



+ - {dh - eli'in)-'^ {Tt{F a F) - tr(^ A ^))} ] , (6) 



where Y^^ — I'^Vy ^ Jy ^ '7"': is the volume of Y, h is the four-dimensional 
part of the 3-form H and 9{x) is a Lagrange multiplier for the anomaly cancellation 
condition, i.e., the Bianchi identity for h. 

From the variation with respect to /i, we obtain dO = 2TT *h. By eliminating h by 
this relation, 9 is promoted to a dynamical pseudo scalar field with the action 



/ *dai A duj *da Ada 

Jxi L 2 2 

A 

J a 



^ a{Tv{F AF)-tr[M AM)} , (7) 
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where fa is the axion decay constant defined by 

f = ^ = = ""pi (8) 

with Vy — , Ui and a are dimensionful axion fields defined by a,j = faCti £^nd a — faO, 
and A is a dimensionless constant 

Here, note that real axion scales for are in general smaller than fa because F*-' ~ 
(is/L)* generally. 

Thus, we obtain two types of scalar fields, the so-called model-dependent axions Oi 
from the internal 2-cycles and the so-called model- independent axion a from B^i,. It 
is clear that the action ([t]) has shift symmetry for both types of fields. Further, the 
model-independent axion a has pseudo-scalar couplings to the gauge and gravitational 
Chern-Simons terms like the QCD axion. In contrast, the model-dependent axions have 
apparently no such coupling and further, they are naively CP-even scalar because it 
is natural to regard the B field as a CP-even field in the heterotic theory. However, if 
we consider the quantum corrections, the Green-Schwartz counter term 

BAXs{F,.^) (10) 

produces pseudo-scalar couplings of the model-dependent axions to Tr(_F A F) and 

B produces a model-independent axion a in type IIB theory and model-dependent 
axions in type HA theory as well, while fields in IIB and a field in HA become CP 
even. In contrast to the heterotic case, these fields do not couple to F A or ^ A 
However, in type II theories, we have various RR form fields Cp, among which C3 in 
HA and C2q{q = 0, 1,2) in IIB theory produce model-dependent axions in the same 
way as those in the heterotic theory. Further, these may be coupled to gauge fields 
via the Chern-Simons coupling of the RR fields with gauge fields on D-branes[lll[^ 

The very important feature of these model-dependent axions is that they are as 
abundant as the non-trial internal cycles. Because the number of such cycles can be 
very huge in the flux compactification of the type IIB theory leading to the landscape, 
a huge number of different axions are expected to be produced in such models [13]. 
Although there have been proposed no definite argument on the Betti number h2{Y) 
in the heterotic models, systematic searches of CY 3-folds in the toric framework 
indicate that 62 (i^) is vary large for a generic CY.[T3] 



2.1.2 Mass spectrum 

If the shift symmetry is not violated at the tree level by flux, branes and compactifica- 
tion (i.e., by moduli stabilization), it can be preserved by perturbative quantum correc- 
tions (for supersymmetric states) . Then, axions acquire mass only by non-perturbative 
effects (possibly associated with SUSY breaking), such as instanton effects as in the 

^It appears that these axions do not have a coupUng to Sf, at least in the lowest order in a' . 
This feature might be used to probe the background string theory of our universe through axions. 
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Figure 1: Various axiverse windows. 

case of the QCD axion. If a light QCD axion really exists, it is natural that there 
survive lots of other light axions coining from the large number of non-trivial cycles in 
extra-dimensions as discussed above. 

Now, let us give a rough estimate of the mass of such axions. In general, the action 
of an axion whose mass is generated by instant effects can be written 

^^-IfaidOf-^'Uie); A^^mV^, (11) 

where 5* is the instanton action. From the relations 

ml,^g-'L%\ fl^9-^LH-\L^)-^=g-^LH-\ S^l-^L^ (12) 

it follows that fa = m-p\/ S. Hence, we have 

ma - A^fa - (MV™pi)5e-^/2 (13) 

For the QCD axion, the total potential is the sum of the QCD contribution and 
the stringy contribution given above: 

Requiring this stringy effect to be less than the QCD instanton effect leads to the 
constraint 

l\,(-4 -S 

a « —^r-F^ < 10"^° ^ 5 w 200 ^ /a w lOi'^GeV, m„ - lO^^^eV (15) 

Thus, it is expected that there are lots of superlight axions whose mass spectrum is 
homogeneous in log to, producing the axiverse. 

2.2 Axion cosmophysics 

As we saw above, typical masses of stringy axions are expected to be much smaller than 
the standard QCD axions. This smallness of the axion mass fi makes the Compton 
wavelength l//i, i.e., the minimum coherence length, comparable to cosmophysical 
scales and opens up possibilities for string axions to provoke the following interesting 
cosmophysical phenomena: [1] 
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Figure 2: Density parameter of axions as a function of /i and fa 



• Birefringence of CMB polarization 

• Step structures in the cosmological power spectrum 

• Black hole instability/bosenova 

• Circular polarization of primordial GWs 

• Penetration of the GZK type barrier of CMB for high energy gamma rays. 

For example, a coherent axionic field acts as a quintessence-type dark energy when 
H <C 3H where H is the cosmic expansion rate, while it behaves as cold dark matter 
when H becomes smaller than ^/3. Typical mass scales are 

fi = 3H = 4:.5x lO-^^eY{H/Ho) (16) 

where Hq is the present Hubble constant. The masses corresponding to 3H at the 
CMB last scattering and the radiation-matter equipartition time are around 10~^®eV 
and 10"^'' eV, respectively. Hence, if an axionic field with mass ^ ~ 10~^®eV comprise 
a non-negligible fraction of dark matter, it produces a step-function-type deformation 
of the CDM perturbation power spectrum [TJ [TS], which may be observable by future 

experiments. When an axionic field with 5 x 10~'^'^eV ^ /i ^ 10~^^eV is a non- 
negligible component, it produces rotations of the polarization of CMB [16] which can 
be observed by the on-going and future B-mode CMB experiments. 

Here, note that for string axions, there exists no universal relation between the 
axion mass /i and the axion decay constant unlike the QCD axion. Hence, we 
can treat them as independent parameters. Then, if we neglect corrections on /i 
by interactions with matter such as the temperature correction, the present density 
parameter fta is less than the DM bound Hum ^ 0.25 for a wide parameter region 
as shown in Fig. [2] The axion abundance in this figure is calculated under the 
assumption that the axion field a has an amplitude of the order of fa initially, i.e., 
during inflation, and no dilution occurs after reheating. Hence, if we assume that the 
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Figure 3: Black hole spacetime 



initial value of a in our observed region is much smaller than that typical value for the 
anthropic reason or else, the present abundance can be much smaller than the value in 
the figure. For example, if fa is around lO^^GeV, as suggested by the argument in the 
previous subsection, the mass range 10~^"eV ~ ~ 10~^^eV is allowed only under 
this anthropic assumption. For /i 3> 10~^^eV for which 3i? becomes smaller than 
before BBN, dilution by decay of heavy moduli may reduce the present abundance of 
the axion. 

Another important mass range comes from the horizon size of astrophysical black 
holes. As we see later, an axion field around a rotating black hole becomes unstable if 
its Compton wavelength is comparable to the black hole horizon size: 



Because the mass M of astrophysical black holes is in the range 

axions with mass in the range 10^^"eV ^ ii ^ 10^^"eV can really produce such 
instability around black holes in binary systems and at galactic centers. In the following 
sections, we study this problem in more details. 

3 Black Hole Basics 

In this section, we briefly overview the basic concepts on black holes that are relevant 
to the superradiance instability. 

3.1 Definition of a black hole 

Let ^ be a weakly asymptotically simple spacetime and be its conformal infinity [17] Din 
order to avoid the appearance of singularity outside a black hole, we require that the 
spacetime is asymptotically predictable from a Cauchy surface E, i.e. .y C D(T,) in 
Under this condition, we define a horizon as the boundary of the region that can be 
observed by the infinity as 




(17) 



H+ 



d{j-[.y))nj+{j). 



(18) 
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Figure 4: Killing horizon 



Figure 5: Horizon of a rotating 
black hole 



Then, the black hole region is defined as the region that cannot be seen from the infinity 
as 

^ = ^ - (19) 

and the region outside the horizon is called the DOC (Domain of outer communication): 

DOC = j-(y,.^)nj+{y,^). (20) 

3.2 Killing horizon 

3.2.1 Stationary spacetime 

A spacetime ^ is said to be stationary if there is a Killing vector ^ that is timelike in 
some region. The metric of a stationary spacetime can be written 

ds^ = -e^^^^\dt + A{x)f + g,j{x)dx'dx^, (21) 

where x = (x*) is the spatial coordinates. The Killing vector ^ can be written ^ ~ dt 
in this coordinate system, hence the corresponding 1-form is given by 

^^=-e'^'^{dt + A{x)). (22) 

The rotation of the Killing vector is defined as 

*(C* AdC*) = -e^^ *„dA. (23) 

A spacetime ^ is said to be axisymmetric if there is a Killing vector field 77 whose 
orbits are all closed and are spacelike in some region. In this article, we mainly consider 
a stationary and axisymmetric spacetime. From the rigidity theorem for black holes, 
a stationary black hole is always axisymmetric if the spacetime is analytic. 

3.2.2 Killing horizon 

A null hypersurface Jif in a stationary spacetime is called a Killing horizon when there 
is a Killing vector that is parallel to the null geodesic generators on A horizon of 
an asymptotically simple and static spacetime with respect to infinity ^ is a Killing 
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Figure 6: Penrose diagram of the Schwarzschild black hole 

horizon if the spacetime is asymptotically predictable and the time translation Killing 
vector ^ is timelike in a neighborhood of J^. 

The black hole of a stationary spacetime is said to be rotating if the time-translation 
Killing vector is spacelike on the horizon. From the rigidity theorem, the rotating black 
hole horizon is a Killing horizon. 

For a stationary and axisymmetric spacetime with a Killing horizon J^, let ^ and 
rj be the corresponding Killing vectors. Then, a tangent vector of the null generator 
of Jif can be uniquely written as 

k = C + nhT]- (24) 
flh is called the angular velocity of the horizon. Further, on Jif, we have 

Vfcfc = Kk <^ Vk^ = -2Kk (25) 
The coefficient k is called the surface gravity of the black hole. 

3.3 Examples 

3.3.1 Static black hole 

A stationary spacetime ^ with the time translation Killing vector ^ is said static when 
the rotation of f vanishes. When a spacetime (^, g) is static, from the rotation free 
condition, we can find a coordinate system locally in which the metric can be written 

ds^ = _e2t^(^)di2 g^^{x)dx'dx^ . (26) 

The simplest and most important example of a static black hole is the spherically 
symmetric black hole solution to the vacuum Einstein equation, whose metric is given 
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by 

ds' = -f{r)df + + rHal; /(r) = 1 - (^)""' - Xr^ (27) 

where dcr^ is a metric of the 71-dimensional unit sphere S"", and A is a constant related 
to the cosmological constant by A = n(ji+i) ' asymptotically flat vacuum case 

with A = 0, this is the unique regular static black hole solution in four and higher 
dimensions from the uniqueness theorem. 

For this spherically symmetric spacetime, the horizon is a Killing horizon and its 
location is given hy r — rh in terms of a solution to f{rh) — 0. The horizon is obviously 
homeomorphic to K x 5". 



3.3.2 Kerr black hole 

In the asymptotically flat vacuum case in four dimensions, a regular rotating stationary 
black hole solution is unique and given by the Kerr solution with the metric 

ds^ = -^de + ^^^^{d4> - ndtf + p^{^ + de^^ , (28) 

where 

A = r2-2Mr + a^ = + cos^ 6*, T = {r^ + a^^ - a'^ Asin^ 9, n = 



2aMr 



r 

(29) 

The global horizon of this spacetime is again a Killing horizon, and its location is 
given by the largest solution to A{r) = as r = rh = r+ = M + (Af^ — a^)^/^. Note 
that A(r) = is equivalent to the condition that the Killing orbit spanned by dt and 
becomes null. The horizon is topologically M x S"^. Because the rotation of the 
time translation Killing vector dt does not vanish, 

uj^di j^O, (30) 

this black hole is rotating. Hence, the time-translation Killing vector is spacelike on 
the horizon as in FigjSj and there appears a region called the ergo region where gu > 
as is seen from p^gtt — o? sin^ 9 ~ A. The existence of this ergo region plays a crucial 
role in the superradiance instability of a rotating black hole discussed in the next 
section. The angular velocity of the horizon Jl/j is determined by the condition that 
dt + Vlhd^ is a null vector as 

2aMrh a a 



4 Bound States and Scattering 

4.1 Particles around a black hole 
4.1.1 Schvifarzschild black hole 

The behavior of a particle around a black hole depends both on the rotation of the 
black hole and whether the particle is massive or massless. 
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Figure 7: The effective potential for a particle with L = 0, • • • , 5 around the 4D 
Schwarzschild BH 



For example, geodesies around a Schwarzschild black hole with the metric (27) can 
be determined by solving the first-order ODE system obtained from the energy and 
angular-momentum conservation laws, 

E = —u ■ ^ — —ut — f{r)i, L = u ■ rj = — r^cf), (32a) 

~e = -fe + j+r^\ (32b) 

where e = 1 for a massive particle and e = for a massless particle. In particular, their 
qualitative behaviour can be easily found from the behaviour of the effective potential 
V{r) in the effective energy conservation law expressed in terms of the r coordinate 
(see Fig.[7|: 

f^ + V{r)^E^- V(r)^(e+^\f{T). (33) 



In particular, we see that there exist stable bound orbits for a massive particle while 
there exists no such orbit for a massless particle. Note that in five or higher dimensions, 
even a massive particle has no stable bound orbit. 

4.1.2 Kerr black hole 

For a Kerr black hole, we have to use the Carter constant in addition to the energy 
and angular momentum to reduce the geodesic equations to a first-order system of 
ODEs in general. However, orbits on the equatorial plane can be determined by the 
energy and momentum conservations laws 

E = -gut - 5t00i L = g^ti + 5000, (34a) 

-e = 5tt^^ + 25t^i0 + 5000^ (34b) 

or the effective potential for the r coordinate, 

r^ + V{r)^E^- V{r)^^ . (35) 
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Figure 8: The effective potential for particles with L — 3 (corotating) and L — —3 
(counter-rotating) around the Kerr BH with a — 0.999. r* is the tortoise coordinate 
defined by dr* = (r^ + a^)dr/A. 



where the particle is corotating with the black hole for L > and counter-rotating for 
L < 0. As we see from Fig. [8] the behavior of a particle is largely different for the 
corotating case and for the counter-rotating case. In particular, the centrifugal force 
on counter-rotating particles is weaker than that for corotating ones. 

4.2 Massless fields around a Kerr black hole 
4.2.1 Flux conservation law 

As we saw in the previous subsection, a massless particle incident to a black hole 
is simply scattered off or absorbed by the black hole, and nothing peculiar happens. 
However, if we consider the scattering of a massless scalar field, an interesting new 
phenomenon happens, when the black hole is rotating. This is because a wave behaves 
as an ensemble of particles. 

To see this, let us consider the free Klein-Gordon field around a stationary and 
axisymmetric rotating black hole satisfying 

D^'D^(b = 0; D^^d^- iqA^, (36) 

where we have considered a charged black hole with the electromagnetic potential 
and a charged field with charge q for generality. Then, we find that the Klein-Gordon 
inner product 

N{<j>,,cl>2) [ (0ii?^02 - (5^0i)</'2) (37) 

does not depend on the choice of a Cauchy surface S in the DOC. In particular, if 
we consider the scattering problem as shown in Fig|9j we obtain the following relation 
among the flux coming from infinity the fiux absorbed by the black hole I jp+ 

and the flux escaping to infinity under the assumption that no fiux is coming from 

Ijf- = 1^+ + (38) 



12 



Figure 9: Scattering of an incidental wave by an AF black hole 



4.2.2 Superradiance 

Let us evaluate each of these fluxes. First, because the wave behaves at infinity as 

(/) « f du;J2- (^"e-*""- + ^+6+*"^"+) e™^; u± = i T / dr/f, (39) 
the fluxes from and to infinity can be expressed as 

=t fdu± f dn^ lim r2(09„^0) = J2f dLu^{\At,J^)s-, (40) 
J Js^ ,„ J 

where {Q)s2 is the average of a function Q on S*^. 

Next, near the future horizon .j/zf^, the wave behaves as 

= 0(r^5/)e»".'-*e-»"*"++""*' w C(6l)e-*"-''++™^, (41) 

where :— lu — mflh, (p — ip — flht, and v-^- = t + j dr(r^ + a^)/A. Because, w+ 
and (f are regular coordinates around this implies that C should be a bounded 

function of 9. Hence, the flux crossing can be calculated as 



= ^J dv+ j^^ d'^-^a l^(d,^ + 2iq^)(^ 
= j du:{uj,-q<^h){rl + a'){\C^,^\^)s-. (42) 



Inserting these expressions for fluxes into (38), we obtain 



^{\AZ,n?) = ^(|A+„J2) + {uj- mn^ - 9$,,) {rl + a'){\C^^rn\') ■ (43) 

If we define the transmission rate T and the reflection rate R of the wave by the black 
hole as 

T := /^+ //^- , i? 11^- , (44) 
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Figure 10: The Penrose process in the ergo region 



we obtain R > 1 because T + R = 1 and T < when the condition 

- q^h = w - mfi/j - q^h < 



(45) 



is satisfied, where is the electric potential of the black hole. That is, the scattered 
wave is amplified and have a larger flux than the incoming wave. This phenomenon is 



called the superradiance and the condition ( 45 1 is called the superradiance condition 



4.2.3 Penrose process 

This curious phenomenon is closely related to the Penrose process in the ergo 
region. In this region, because the time-translation Killing vector ^ is spacelike, the 
energy E = —p ■ ^ with respect to the spatial infinity can become negative even for 
a future-directed timelike 4-momentum p^. Hence, if a particle incidental to the ergo 
region decays into two particles there, the outgoing particle can have a large energy 
than the incidental one if the other decay product has a negative energy as shown 



in Fig 10 This is called the Penrose process. The produced negative energy particle 
cannot go outside of the ergo region, and is eventually absorbed by the black hole 
to reduce its mass. Hence, the Penrose process makes it possible to extract energy 
(and angular momentum) from a rotating black hole. The superradiance above can 
be regarded as being produced by the same mechanism because the superradiance 
condition can be written 



k ■p> 0, Pf_,(f> = i-idi_, - qA^ 



(46) 



where k = dt + ^hdip is the future-directed null generator of the horizon This 
equation implies that p is a past-directed timelike vector because the physical momen- 
tum p is timelike. 
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fiM = 0.5, / = 1 



fiM = 0.9, I = 2 



Figure 11: The effective potential for a massive scalar field around the Kerr BH with 
a = 0.999. 



5 Superradiance Instability 

Superradiation does not cause any physical problem by itself. However, if we consider 
the gedanken experiment to put a rotating black hole inside a box with reflective 
boundary, superradiance provokes an instability called a black hole bomb [El [2Qj |2T] 
because amplified waves of a massless field by superradiance are reflected back toward 
the black hole by the surrounding mirror wall of the box, and the repetition of this 
process produces an unbounded growth of the waves. This instability lasts until the 
central black hole loses its whole angular momentum. 

In reality, of course, it is impossible to put a black hole in a box, but Damour, 
Deruelle and Ruffini[22 pointed out that we can realize such a situation practically 
just by considering a massive field instead of a massless field around a rotating black 
hole. This is because a massive field can have bound states and the effective potential 



provide an effective outer wall for the field as shown in Fig 11 In this section, we look 
at the instability of such a massive scalar around a Kerr black hole in details. 



5.1 Massive scalar equation around a Kerr black hole 

The field equation for a massive free scalar field 

(□ - = 

can be separated by 

$ = Rim{r)Sim{0) eyip{-iujt + imcj)), 
into the ODE for the angular mode function Sim{0), 
1 



sin6'd6' 



d ( . „dSi 
sm 



d9 



a^{uj^ - cos 



sin^ e 



A 



Im 



Sim — 0, 



(47) 



(48) 



(49) 
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with the separation constant A/„i, and the ODE for the radial mode function Ri„i{r), 



d f ^ dRim 
dr \ dr 



-{lu^o^ + fi^r^ + Aim) 



Rim — 0. 



(50) 



The equation ( 49 ) for the angular mode function, as an ODE with respect to the 



independent variable x — cos 9, have two regular singularities at a; = ±1 with index 
±m/2 and an irregular singularity at infinity. Hence, it has a series of regular solutions 
Sim = 5'["(a;;c) with c = a(w^ — fJ-^Y^"^ on the interval —1 < x < 1 corresponding to 
discrete values of Aim labeled by an integer ^ = 0, 1, 2, • • • . In the limit c 0, it 
reduces to the standard associated Legendre function as 



sr^pr, Aim ^1(1 + 1) 



(51) 



The equation ( 50 1 for the radial mode function can be rewritten in terms of w = 

(52) 



(r^ + a^)^/^Rim as 

d^U r 9 , >n 

which defines the effective potential 

fi'^A AamujMr ~ a^rr? + A[A;„ + [b? - [p-)a^ 



V 



+ (r^ + a^)^ 

A{2Mr^ + a^r^ - AMa'^r + a^) 
^ (r2 +a2)4 ■ 



(53) 



This potential always approaches /i^ at r = oo and lo^ — lu'^ at horizon as illustrated in 
Fig. |11| In order to study the superradiance instability, we have to look for a unstable 
solution to this equation satisfying the following boundary condition: 



At infinity 
At horizon 



Rlr. 
Rl.r. 



Im +ikr* 



k 



2U/2 



(54a) 



Cime"''^''' 6"*'^*+™'^. (54b) 



Note that the first condition implies that no wave is coming from infinity and the 
second condition requires that waves are purely infalling at the horizon. 

5.2 Instability condition 



From the field equation, for the solution of the form (48), we obtain the following 
energy integral: 



= i 



J desm0{~uj*)A<i>*dr^ 



dr / dO sin t 



p'\u;\\u:g''-2mg'*m 



+uj*p' {jn'g^^ + |$|2 + ^* (A|9,$p + 



(55) 
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If the solution satisfies the boundary condition (54a) and (54b), this identity can be 
written in terms of cj = ujji + iuji and R — Rim exp(iw,tr*) as 

BujR{mQh - i^Rj + t^/(Ci + C2i^i) = Auoi, (56) 



where 



B^{rl+a')^-'^-\R\l^^^, 
Ci = y dr2re2"^"+|J?|^ 



02=0^ 



dre2"^"+ j d9sm^9\RSp\^, 



and 



A = I dre^'^"'+ I d9sm9 
sin 9drR 



\R\^\deSl^ 



o'i-gtt)\drmsn' 



l^n — TT-. — T^R 



sm 9 



2 2 



ii?5ril37) 



where Q is a real quantity and P is 



P = p\T'^a'f{p'^^a'^\^^9){guf 

a^Mr^ sin^ 9 



+8Ma^ sin^ 9 



[p^ + 2a^ sin^ 6*) . (58) 



Obviously, A,B,Ci and C2 are all positive and finite if uj'^ < p? , which guarantees 
that the wavefunction falls off exponentially at infinity, and if the contribution from 
the integral in the ergo region where {—gtt) < is not dominant. In that case, from 



the identity (56), it follows that w/ > when the superradiance condition is satisfied. 



This implies that the scalar field grows exponentially in time as $ oc exp(a;7t). If any 
of these conditions is not satisfied, we cannot draw a definite conclusion about the 
instability. Thus, we obtain the following sufficient conditions for the occurrence of 
superradiant instability [23]: 

i) The mode is bounded. 

ii) The wavefunction is peaked far outside the ergo region. 

iii) Lo is nearly real: ^ cj^. 

iv) Lo satisfies the superradiance condition: w/? < m^lh- 



5.3 Growth rate 

In this section, we estimate the growth rate of the superradiance instability of a massive 
scalar field with mass p around a Kerr black hole with mass M by three different 
methods. The first is based on the WKB approximation that is valid when Mp ^ 1 
in the absolute units c = h = G = 1. The second is the matched asymptotic expansion 
method that is valid when Mp ^ 1. Unlike these semi-analytic methods, the third is 
a purely numerical method based on the continued fraction that was first introduced 
by Leaver [23]. In this section, we adopt the absolute units unless otherwise stated. 
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Figure 12: Division into four regions for the WKB approximation 



5.3.1 Large mass case 

In the absolute units, the typical scale of the background geometry is M, while the 
Compton wavelength of a particle with mass ^ is l//i. Hence, when M/i 3> 1, the 
wave length of the scalar field is much shorter than the geometrical scale, and the 
WKB solution provides a good approximation for the wavefunction. The growth rate 
of the superradiance instability with this approximation was first estimated by T.J.M. 
Zouros and D.M. Eardley[23]. We follow their arguments but some numerical error in 
their result is corrected. 

Because we are interested in the superradiant instability, we consider the quasi- 
bound state with < /i^. Because the instability rate becomes larger as the wave 
amplitude at horizon is larger, we only consider the case a; a; ~ /i in which the tunneling 
transition late through the potential barrier becomes maximum. For such a value of 
lor, the range of the coordinate r is divided into four regions I(r < ri), II(ri < r < r2) 
III(r2 < r < r^) and IV(r > rs) according to the sign of — V{r) as shown in Fig 12 



In the oscillatory regions I and III where > V{r)^ the WKB solution for the radial 
mode function u = (r^ + a^Y^^Rim can be written in the form 

where A± are constants and 

e(r) = / k{u)du, k{r*) = {lu^ - V{r)y^'^. (60) 



Here, Tq is the r* coordinate of a reflection point where — V{r). Meanwhile, in the 
regions II and IV where cj^ < V{r), the WKB solution reads 

u = K(r*)-i/2 {B_e-^('") + B+e^('')} , (61) 
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where B± are constants and 



/(r) 



r K{u)du, K{r*) = {V{r) - u;2)i/2_ 



(62) 



Because these WKB solutions diverge at the reflection points where k and k vanish, 

we cannot directly connect these WKB solutions in different regions by the standard 
regularity requirement. The standard method to cope with this situation is to utilize 
the exact solution around each reflection in the case where the potential is locally ap- 
proximated by a linear fimction of r*, which can be written in terms of Airy functions. 
In the oscillatory region, it can be written in terms of the Bessel functions J±i/3 as 



u 



1 



'2-Kk 



{c+Ji/3(|e|) + c7_j_i/3(|e|)} 



and in the tunneling region, as 



u = 



{-C+/l/3(|/|) + C_/_V3(|/|)} 



1 



(C_ - C+)el^l + (C_e-* - C+e-"^) e'W 



(64) 



By applying this method to connecting adjacent regions starting from the region 
1 with the inf ailing boundary condition A^_^_ = 0, we obtain the following relations in 
order: 



igii ^ e-'^'/'^Ao, B^} = 0, 

5IV ^ g-3^V4e/n+ieni^^j^ 



(65a) 
(65b) 
(65c) 
(65d) 



and 



where 



(e^'/3 _ i)_Biv ^ 2e5^'/i2e-f" cosem. 



K{r)dr*, 6111= / k{r)dr* 



(66) 



(67) 



Here, we have taken the base point for each phase integral / or at the left end of each 
region except for the region I for which r = ri is taken. Because we are considering 
a bound mode with cj^ < /i^, B^J^ should vanish. This leads to the Bohr-Sommerfeld 
quantization condition for the bound state frequency, 



k{r)dr* 



n + - ] TT, 



n = 0, 1,- 



(68) 
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Figure 13: Flux integral to estimate the growth rate 



Now, we can estimate the instability growth rate using this WKB approximation 
for a quasi-bound state solution. We first note that the KG inner product for an exact 
complex solution (cx exp(— zwt)) representing an unstable mode has to vanish if it is 
estimated on a hypersurface Eg with the inner boundary at the bifurcating 2-sphere 



(see Fig 13). It is because the boundary of the hypersurface does not move by a time 
translation and as a consequence, the KG inner product must be time- independent, 
while it has to be proportional to exp(2cj/f). This implies that the integrand is negative 
in a region close to the horizon on such a hypersurface. This can be confirmed by look 
at the explicit expression for the KG inner product, 

iV($,$) = e2-^*y {ujR{{r^ + a^f-a^A(:)~2maMr}\Ri^\^"^, (69) 

where C is the average (sin^0|5™p) over 5^. 

Thus, this type of a hypersurface cannot be used to estimate the instability growth 
rate. We therefore take a hypersurface E that has the inner boundary at the future 
horizon and is parallel to the infalling null direction inside the potential barrier as 
shown in Fig. |13| Then, because the wave is almost purely infalling inside the potential 
barrier due to the boundary condition at the horizon, the flux integral on this part of 
the hypersurface E can be neglected. Further, we can confirm that the peak position 
of the barrier is around the boundary of the region where the flux in the Vi direction 
is negative. Hence, by translating this surface by time translation, we obtain from the 
flux conservation 

-u4rl+a^)\R\^ ^2iOiN^,i^,<^). (70) 
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From this, w/ can be estimated as 



-76 



-2/,, 



7 



dr* 



4cos^ 



2maMr 



(71) 



(72) 



Zouros and Eardley numerically estimated the instability growth rate for wide 
ranges of the parameters a/M, iJ,M,l,m and u! using this formula. They found that 
the growth rate becomes the largest for i) the smallest I, ii) the largest possible m, 
i.e., m = I, iii) the largest possible a/M, i.e., a/M ~ 1, and iv) the largest possible 
uJn, i.e., ujji ^ 0.98/i < mVLii, and gave the estimate for the maximum growth rage for 
a given /xM, 

Muji - 10"'^ exp {-l.M^iM) , (73) 



The very tiny prefactor 10 ^ comes from 7 in (71). 



5.3.2 Small mass case 

In the limit /iAf <C 1, we can estimate the instability growth rate by a different 
method[25j(Cf. |26j). This method utilizes the fact that we can approximate the 
mode function Rim by known analytic functions in two asymptotic regions. 

First, in the asymptotic region r ^ M, we can approximate the ODE (50 1 for 
R = Rim{r) by 



(P{rR) 



2Af^2 ;(/ + !) 



(rR) « 0, 



(74) 



which has exactly the same form as that of the Schrodinger equation for a hydrogen 
atom. Hence, when = ~lo^ > 0, it has a sequence of quasi-bound state solutions, 



R=-W. 



/2{x) 



{x = 2ar > 1), 



z/ = M/x^/(T = / + n+l + (5t/, (n = 0, 1,2,---), 



(75) 
(76) 



where 8v is a complex number representing the deviation from the exactly hydrogen- 
type wavefunction. This solution behaves in the region aM <C a; ^ 1 as 



i-i 



Next, in the region /ir <C I, the ODE (50 1 can be approximated by 
d 



z{z + l) 



dz 



z{z + l) 



dR 

dz 



+ {P^ - l{l + l)z{z + 1)} R ^ 0, 



where 



r — ra 



P = 



2Mr. 



+ 



This can be exactly solved, and the solution infalling into horizon is given by 

iP 



R^C 



z 



z - 



1 



F{-l,l + 1,1 + 2iP]~z). 



(77) 

(78) 
(79) 

(80) 
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Figure 14: The instability growth rate for I — m ~ 1. 



This solution behaves in the overlapping region 1 <C z <C 1/{ujrM) as 

^ r20!r(i + 2*p[ , ,+i_nr(i + 2zF^ , 1 
urii + i + 2iP)^^ ' {2i + iy.r{-i + 2iP) ' 

Matching the two approximate solutions in the overlapping region, we obtain 

1 2 ' 



(81) 



Si' = 2iP [2a{r+ - 



i+i i2l + l + n)l 



n 



{2iy.{2i + i) 



This determines the instability growth rate as 



1/2 



where 



7 = 



2*^+^21 + 1 + n)l 
n!(; + l + n)2'+4 V(20!(2^ + 1) 



II 



Uif+iP^). (82) 



(83a) 
(83b) 



Clearly, for I = m = 1 and a/7\f ^ 1, the growth rate takes the maximum value 

a 



24Af2 



{piMf. 



(84) 
(85) 
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5.3.3 Numerical estimation 



The estimations of the instabihty growth rate, (73) for /iAf ^ 1 and (851 for jiM <C 1, 
suggest that it becomes maximum at around fiM ^1. To prove this and find the 
maximum growth rate, we have to solve the ODEs for mode functions numerically. 
This numerical study was first done by Cardoso and Yoshida[17] using the continued 
fraction method [53] and later extended to a larger parameter region by DolanPSl. This 
method has been also used in calculating the quasi-normal mode frequencies of higher- 
dimensional black holes[27| |29] on the basis of master equations for perturbations [30l 
[211 [3H |33J 133] and establishing the SR instability of simply rotating adS black holes[35l 
136] . We briefly overview this numerical method in this subsection. 
He first expanded the radial mode function R as 



Rir)^ . " w-i ^"''^E«"^"' (86) 

^ ' n=0 

where 

2r+{u:-mVlh) , 2 2n1/2 (2 o 2n / ^o^n 

q = , (7 = [ji -uj ) ' , X = -[ji - 2uj ) a. 87 



Inserting this into ( 50 ) , we obtain the three term recurrence relation for the expansion 
coefficients a„: 

a„a„+i H- /3„a„ -I- 7„a„_i = 0. (88) 

with 

a„ = (n + l)(n + Co), /3„ = -2n^ + (ci + 2)n -t- C3, 7„ = + (c2 - 3)n + C4, 

where ci, • • • ,04 are constants dependent on a;, cr, to and Kim- 

The point is that under the assumption that a„+i/a„ tend to zero as n — > cx3, we 
can solve this recurrence relation in terms of a continued fraction as 

a„+l 7,1+1 7„+l an+lln+2 a,i+27ji+3 
^I^TTi" - Q a • y^^l 



+3- 



Because ai/oo = —po/ao, this equation with n ^ gives the eigenvalue equation for 
Lo = uji + iuji: 

/3o-f^f^f^--- = 0. (90) 

Pi— P2— P3 — 

The convergence of this continued fraction is very rapid, and by truncating it at 
some appropriate level, we can easily obtain an approximate algebraic equation that 
enables us to determine with good accuracy. As an example, we give the result of 
our numerical calculation for I = m = 1 in Fig|14[ 



6 Axionic Instability of Astrophysical Black Holes 

Now, let us discuss what kind of astrophysical phenomena the superradiance insta- 
bility provokes in realistic astrophysical systems with black holes. Description in this 
section are largely based on the axiverse paper and the paper by Arvanitaki and 
Dubovskv[37j. but it also contains new results obtained by us. 
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Figure 15: The instability strip in the fi — M plane where the instability time scale is 
shorter than the cosmic age. 



6.1 Instability strip 

To start with, let us summarize the estimations obtained in the previous section. First, 
in asymptotic regions, the growth time scale of the superradiance instability of a scalar 
field with mass around a Kerr black hole with mass M and angular momentum aM 
is approximately given by 

r r 10V«4"« ;a<,»l, a = l 

GM 1 24(^)-'(a3)-^ ;a,«l, ^ ^ 

where 



a„:^GMa^ ^ ^„ ■ — . (92) 

" ^ 1.34 • 10-iOeV Mo ^ ' 

It becomes maximum at ~ 1: 

« 0.2 • lO^GM; ag ~ 0.44, a/M ~ 0.999. (93) 

This result implies that if there exists a scalar field with mass /i ^ 10~^°eV in 
nature, a cloud of the scalar field with a large amplitude will show up within one hour 
or so around any black hole with the solar mass due to the superradiant instability 
starting from tiny quantum fluctuations. This implies that the superradiant instability 
may provoke interesting astrophysical phenomena in reality, because the QCD axion 
a, whose existence is high probable taking account of the strong CP problem, can have 
a mass as small as 10~^°eV if the PQ symmetry breaking scale is as large as lO^'^^GeV. 
This mass range is allowed if the initial amplitude during inflation is accidentally much 
smaller than a typical value. 

Of course, this instability does not grow unboundedly because unstable modes have 
corotating angular momenta due to the superradiance condition and as a consequence, 
the central black hole loses its angular momentum by superradiance. Because the 
growth rate of the instability decreases rapidly with the decrease of the black hole 
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Figure 16: The wave amplitude of the most unstable mode(^ = m = 1) on the equa- 
torial plane (the left panel) and on the vertical plane in the most extended direction 
(the right panel). 



angular momentum as shown in FigfM} and the instability stops when the black hole 
loses a non-negligible fraction of its angular momentum. 

The most characteristic feature of this superradiance instability is the sensitivity 
of the growth rate on the masses of the scalar field and the black hole. Due to this, 
if we plot the region in the fji — M plane where the instability growth time scale is 
shorter than the cosmic age (~ 14Gyr), we obtain a narrow strip as shown in Fig. 15 
For example, if there exists a scalar field with mass /i « 10~"'^''eV, only systems with 
a black hole with mass in the range lO^M© — IO^Mq are affected and have smaller 
angular momenta than the other systems. 

6.2 G-atom 

Another peculiar feature of the superradiance instability of a black hole is the quantum 
nature of the phenomena. To see this, first recall that the instability occurs for bound 
states with w « /i. If we approximate its energy levels by the hydrogen-atom type 



ones in the small mass case, from (83a), we have 

„2 



UJ 



R 



1 



where n = n' + I + 1 {n' ^ 0,1,2 ■ 
n' ^ 0, n can be estimated as 



UJR < mnh, (94) 
In general, for the most unstable mode with 



Hence, the mode is peaked outside the ergo region and far from the horizon: 



r 

Rn 



fir 



Aa, 



(95) 



(96) 
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This implies that the most unstable bound states are quantum for near extremal 
cases! Hence, the superradiance instability produces a G-atom (gravitational atom) 
consisting of the central black hole and a surrounding axion cloud in quantum states. 
Fig. [16] illustrate the shape of this G-atom. 



6.3 Gravitational wave emission 

As shown in Fig. |16[ the axion cloud of a G-atom is not spherically symmetric and 
rotating. Hence, it emits gravitational waves. The efficiency of this emission can be 



roughly estimated by the quadrupole formula. First, from (96 1, the cloud radius Tc is 
around Af (/-f l)^/ag, and we can assume that it is approximately Kepler rotating 
with the angular frequency Q — (Af/rj^)^/^. Then, from the quadrupole formula, the 
power of gravitational waves emitted from the axion cloud is 

where N is the number of axion quanta in the axion cloud, e — fiN/AI is the ratio of 
the cloud mass to the black hole mass. 

Here note that in general, the quadrupole formula estimates the gravitational wave 
emissions due to a slow change in the trajectories of source objects, which corresponds 
to transitions of the axion cloud from higher energy level to lower ones in the present 
problem. This implies that levels with / > 2 are relevant. Hence, the time scale for 
the axion cloud to lose a fraction e of the gravitational energy of the black hole by 
quadrupole-type gravitational wave emission becomes 

From this, we obtain 

TGW/Tsr « 0.1 X e-l-844(«.,-2) (2/^^)10 (99) 

This estimation implies that for ag < 2, this time scale is much longer than the SR 
instability time scale, hence the quadrupole gravitational wave emissions do not stop 
the growth of the G-atom by instability. In contrast, for ag > 2, tqw becomes shorter 
than the SR instability time scale for e ^ 10^^. This implies that for this parameter 
range, the SR instability stops growing when the cloud mass becomes eM where the 
value of e is determined by the condition tqw = ''"SR because axions in unstable levels 
go to stable levels, provided that the other processes do not affect this balance. 



From (97), the observed amplitude of the GW metric perturbation is estimated as 



— -(T^)(c^)m(ToSg)(?)"(<^)^ 

This indicates that massive black holes at galaxy centers produce gravitational waves 
with intensity observable by the next generation GW detection experiments such as the 
advanced LIGO, if there exists axions with mass in the range lO-i^eV ^ fi^ 10-20eV. 

Of course, the estimation by the quadrupole formula may be too crude, and more 
precise estimations taking account of the quantum nature of the axionic cloud are 
required. Some preliminary analysis has been done in [|37)]. but a systematic study is 
yet to be done. 
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Figure 17: The fate of an axion cloud around a black hole formed by instability. 



6.4 Bose nova 

In the discussions so far, we have neglected axion self-interactions. Because the axionic 
field grows coherently by superradiance instability, this approximation may become 
bad when the amplitude of the axionic field becomes large. 

In order to estimate the effect of self-interactions, let us consider the axionic field 
d) with the action 



S 



(fx. 



1 



(V^)^ 



Wfa) 



(101) 



Because the bound states of a G-atom is non-rclativistic, we assume that the axionic 
field can be approximately written in terms of a slowly varing function -0 as 



1 



Then, when \4>\/ fa ^ 1, we obtain a non-relativistic effective action 



(102) 



^1 



NR. 



2/u 16fi 



(103) 



The last term is the leading term of the interactions for small amplitudes. Because 
the corresponding energy is negative, it corresponds to an attractive force. Hence, it is 
expected that the axion cloud of the G-atom collapses due to this attractive force when 
the cloud becomes dense enough like the bosenova phenomena in the Bose-Einstein 
condensate as discussed by Arvanitaki and Dubovskv|37). 

To see this, we approximate the total energy of the axion cloud as 



E = 



V 
2(1 
N 
2^ 



1 



agN 



V 
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Figure 18: The energy of an axion cloud as a function of the cloud size R 

where r is the distance of the cloud from the black hole center, and R is the scale 
describing the extension of the cloud. By minimizing E w.r.t. r, we get the Kepler 
radius of the axionic cloud: 

This energy becomes maximum w.r.t. R at R — Rm- If Rm < ^c, the cloud with 
R r ^ Tc is stable, while if Tc < Rm, E becomes a monotonically increasing function 
of R in the range R < rc as shown in FigjlSj and the cloud initially with R ^ r becomes 
unstable and collapses rapidly. This happens when the total mass of the axion cloud 
exceeds the critical value given by 

;2 j-2 



v2 }->n 2 



<R^ ^ m7V>^-^ ^ ^ = ^ > « 10"'- (105) 



Here, it should be noted that the real self-interaction of the axion field is not quartic. 
When the amplitude increases, the higher-order terms with alternative signature come 
to contribute, and eventually when |'/>|//a becomes of order unity, the full potential 
in proportion to cos((/)//i) should be used. Hence, it is a very interesting problem to 
find what really happens. To see this fate of a growing G-atom, we are now studying 



the evolution of an axion cloud around a Kerr black hole using the exact action ( 101 ). 
According to the results we have obtained so far, the bosenova collapse of the G- 
atom really occurs when the amplitude of (f) becomes comparable to fa, i.e. when 
e ~ 10~^. After the collapse, the large fraction of the energy deposited in the axion 
cloud falls toward the black hole as waves. Figure [19] show two snapshots of the field 
amplitude on the equatorial plane before and after the bosenova collapse for a 3D 
numerical simulation of the axion evolution starting from the 2p mode for which the 
SR instability becomes maximum. The details of the results and physical analyses of 
them will be published elsewhere. 
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Figure 19: Snapshots of a numerical simulation for bosenova collapse 
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7 Summary and Discussion 



In this article, we have overviewed the basic idea of the axiverse and its cosmophysical 
imphcations in order to show that superhght string axions can provide a new cosmo- 
physical tool to probe the string theory/M-theory as the ultimate theory of nature. In 
particular, we have picked up the superradiance instability of astrophysical black holes 
as one of the most fascinating cosmophysical phenomena provoked superlight axions 
and given detailed accounts of its mechanism, the estimation of the instability growth 
rate and observational consequences. 

Although many of the basic ideas of the axiverse and its cosmophysics are not new, 
the recognition of their importance is rather new. Furthermore, new ideas about axion 
cosmophysics are proposed every week. Therefore, there remain lots of problems to be 
studied systematically. For example, in order to determine what really happens when 
the superradiance instability grows, in addition to gravitational wave emissions and 
the non-linear dynamics of axion clouds, we have to calculate various other processes 
such as the direct axion emissions due to non-linear interactions of axions and the ra- 
dio emissions by the Primakov-type process taking account of the ubiquitous existence 
of strong magnetic fields in active black hole systems. In order to construct a model 
describing a realistic system, we also have to take into account the energy and angular 
momentum supply to the black hole by accretion. Along with these cosmophysical 
investigations, it has also a crucial importance to construct string compactifications 
providing an effective four-dimensional cosmology consistent with all observations and 
experiments and calculate the axion spectrum on the basis of such compactifications. 
We can extract information on the ultimate theory only through the comparison of 
the string theory predictions obtained by such a first principle approach and the cos- 
mophysical observations. Systematic investigations of both of these string theory and 
cosmophysics problems are now challenged under "Axiverse Project" in Japan. 
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